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Abstract 



The skin of Tursiops Truncatus is corrugated witli small, quasi-periodic ridges running 
circumferentially about the torso. These ridges extend into the turbulent boundary 
layer and affect the evolution of coherent structures. The development and evolution 
of coherent structures over a surface is described by the formation and dynamics of 
vortex filaments. The dynamics of these filaments over a fiat, non-ridged surface is 
determined analytically, as well as through numerical simulation, and found to agree 
with the observations of coherent structures in the turbulent boundary layer. The 
calculation of the linearized dynamics of the vortex filament, successful for the dy- 
namics of a filament over a flat surface, is extended and applied to a vortex filament 
propagating over a periodically ridged surface. The surface ridges induce a rapid para- 
metric forcing of the vortex filament, and alter the filament dynamics significantly. A 
consideration of the contribution of vortex filament induced flow to energy transport 
indicates that the behavior of the filament induced by the ridges can directly reduce 
surface drag by up to 8%. The size, shape, and distribution of cutaneous ridges for 
Tursiops Truncatus is found to be optimally configured to affect the filament dynam- 
ics and reduce surface drag for swimming velocities consistent with observation. 
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Chapter 1 
Introduction 



1.1 Dolphins 

A dolphin gliding through water with grace and speed presents both a beautiful per- 
formance of nature and a persistent enigma to fluid dynamicists. Traditional models of 
fluid flow applied to a swimming dolphin suggest this mammal is capable of aberrantly 
large powers of muscular exertion. However, studies of dolphin energy expenditure 
show that a dolphin has a metabolism similar to that of other mammals [||], implying 
that these animals, rather then having unusual muscular capability, may posses some 
mechanism that allows them to slip through water with unusually little effort. Sir 
John Gray was the first to hypothesize that dolphins may posses an ability to impede 
the development of turbulence over the skin surface and hence reduce the surface 
drag of a swimming dolphin ||2|. A method of reducing boundary layer turbulence 
and surface drag would be extremely useful if it could be implemented efficiently, and 
the actualization of such a method has posed a theoretical and experimental challenge 
to researchers. 

1.1.1 Skin pliability and active control response 

Some researchers have attempted to reproduce the dolphin's conjectured drag reduc- 
tion by studying - and attempting to duplicate - the pliable nature of dolphin skin 0. 
It was thought that the surface drag induced by a turbulent boundary layer could be 
reduced if the surface was able to yield to the microscopic disturbances produced by 
turbulent eddies. Unfortunately, this approach was not significantly successful. More 
recently, a large effort has been extended towards exploring the possibility that active 
control of surface shape in response to the forces from the turbulent boundary layer 
may be able to reduce turbulent development. The proposal that this sort of direct 
active response is implemented by dolphins is ostensibly supported by the existence 
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Figure 1.1: A sketch showing the distribution and orientation of cutaneous ridges 
over the body of a dolphin (From Shoemaker and Ridgway 0). 

of systems of nerves and cutaneous muscle that appear capable of actively control- 
ling the skin surface. However, detailed measurement of dolphins' skin reaction to 
microvibratory stimuli indicate a minimum response time on the order of 10~^s [Q, 
whereas the timescale associated with the passage of coherent structures over the skin 
is on the order of 10~^s This indicates that a dolphin's skin is not able to respond 
to the passage of individual coherent structures. 

1.1.2 Cutaneous ridges 

It is only recently that Shoemaker and Ridgway documented a feature of odontocete 
skin that had been virtually unknown previous to their work; namely, the existence 
of cutaneous ridges covering most of the body surface |^. These corrugations, with 
furrows running circumferentially around the whales, are reported to be /i ^ 10 — 
60/im in height when the animal is at rest, with a separation of L ~ 0.4 — 1.7mm, 
and vary slightly in scale over the dolphin's body. The ridges are absent from the 
dolphin's snout and from much of the head - regions associated with laminar flow, 
as well as from the flippers and tail, (Figure fLTj ). The fact that the corrugations 
exist only in regions of steady turbulent flow is a strong indication that these ridges 
are associated with turbulent development. The absence of the corrugations from the 
flippers and tail indicates that the ridges may not be useful in these regions which 
are exposed to large velocity and pressure variations. 

Microscopic studies of sections of skin over the dolphin's body show the capillary 
and muscular structure to be correlated with cutaneous ridge spacing. This anatom- 
ical analysis is consistent with the possibility that a dolphin may actively tune the 
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Figure 1.2: A photo-micrograph of a cross section of dolphin skin, showing the cuta- 
neous ridges and underlying dermal structure (From Shoemaker and Ridgway Q). 

amplitude of surface corrugations to correspond with swimming velocity. Although 
the measured surface ridge height is small for a dolphin at rest, producing a ridge 
height to length ratio of ~ 0.025, a photo-micrograph of a section of dolphin skin, 
(Figure pT^ ), displays a much larger ratio, ~ 0.2. This indicates that a live dolphin 
may be reducing the ridge height to the smaller value while at rest by using cuta- 
neous muscle, while the "relaxed" ridge height ratio, indicative of the available range 
of ridge heights, is much larger. This hypothesis is further supported by reports from 
trainers at Sea World that the surface ridges become more pronounced with greater 
swimming speed 0. 

The existence of regular corrugations with tunable amplitude covering the bodies 
of these efficient swimmers, in regions associated with steady turbulent boundary layer 
flow, presents a challenging new mystery. This challenge is met in this dissertation, 
which presents a theoretical model for and an analysis of the effect of such a wavy 
boundary on turbulent development. 

1.1.3 Fluid flow development over dolphin skin 

Allthough cutaneous ridges were observed and measured by Shoemaker and Ridgway 
in seven species of toothed whales (odontocetes), ranging in body length from two to 
eight meters ||^, no correlation between body size and corrugation scale is visible in 
the data. It may be possible, through future study, to establish a correlation between 
corrugation scale and swimming speed. However, detailed information regarding the 
swimming speeds of these various species is not currently available. Hence, a cor- 
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relation between swimming speed and corrugation size is not presently possible to 
establish. The best observational data currently exists for the bottlenosed dolphin 
{Tursiops Truncatus), and this species is taken as representative for this analysis. 
Adult Tursiops vary in size from two to four meters, mass approximately 200kg, 



and have been observed to be capable of bursts of swimming speed up to 50¥^, though 



constant cruising speed is reported to be roughly U ~ 10-^ 

Since the boundary layer thickness at the location where the transition to turbu- 
lence occurs is small, approximately 



51 ~ 420- = 420 ^ = .15mm (1.1) 



U W, 

hr 



at cruising speed, the circumferential curvature of the dolphin's body is not significant. 
Thus, a laterally flat dolphin may be assumed for the purpose of studying boundary 
layer turbulence. The transition to turbulence occurs for flow over a flat plate at a 
calculated position of 



' s 

km 
hr 



c^6x 10*— = 6 X 10^ ^ = 2.2cm (1.2) 

U 10^ ^ ' 



from the edge of the plate. However, transition to turbulence in flow over a dolphin 
is signiflcantly delayed by the pressure gradient induced by the longitudinal curva- 
ture of the dolphin's head, and is observed via bioluminescent excitation to begin 
approximately 40cm from the snout 0. 

Allthough transition to turbulence in flow over a dolphin's curved head is delayed 
by the favorable pressure gradient, once initiated, turbulent development and evolu- 
tion is assumed to be of the same character as turbulent development over a flat plate 
with no pressure gradient. 
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Chapter 2 

The boundary layer 



A consideration of the effect of surface corrugations on turbulent flow development 
starts with understanding the evolution of the boundary layer over a flat surface. 

The flow of water, with density p ~ 10^ and viscosity u ~ 10~^ — , is described 
by the Navier-Stokes equation, 

dtlt = dtlt +(lt ^ --"^p + vV'^lt (2.1) 

in which = ^ ~ iQ-^n^ jg the kinematic viscosity, dt is the advective derivative, 
and the fluid velocity fleld, it — u^x + Uyy + w^z, is restricted to be incompressible, 

"^•l^^O (2.2) 

Incompressibility may be assumed, provided fluid velocities remain much less than 
the speed of sound in the fluid. All velocities are described with respect to the rest 
frame of the boundary surface. Therefore, the skin of a dolphin swimming at speed U 
through stationary fluid is modeled as a stationary flat plate immersed in fluid with 
background (free stream) velocity of limj^^oo Ux = U m the x (streamwise) direction. 
The plate is oriented such that y is normal to the plate surface, which lies in the 
plane y = 0, with the leading edge of the plate located at x = 0. The z coordinate 
runs along the span of the plate, with arbitrary origin. 



2.1 The laminar boundary layer 

Viscosity plays the primary role in boundary layer evolution as the irrotational fluid 
flrst encounters the leading edge of the flat plate. The fluid at the plate surface 
is attached to the surface via molecular interactions and must satisfy the "no slip" 
boundary condition, 

^"(0; > 0, y = 0) = (2.3) 
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Figure 2.1: The horizontal velocity profile, ^"^jf^ , of the laminar boundary layer. The 
vertical scale is in units of ^/^. The heights of the displacement thickness, 5*, and 
99% thickness, ^gg, are shown. 



The equations of motion, (|2.1| ), for the flow may be solved numerically or approx- 
imately to obtain the laminar boundary layer profile, Ux{x,y), (Figure [2Tl| ). The 
profile may be characterized by several descriptors of boundary layer thickness. The 
99% thickness, ^gg, is the most common descriptor, and is defined arbitrarily as the 
height at which Ur^^y = Sgg) = .99U. Another less arbitrary descriptor, the displace- 
ment thickness, 6*, is defined to be equal to the distance the wall would have to be 
displaced into the fluid such that the mass flux of fluid over the boundary would be 
the same as that for fluid flowing over the wall "free slip" with the free stream velocity 




and may be considered the middle of the boundary layer. 
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The boundary layer fluid induces a shear stress of 



y=0 



on the wall in the x direction. This stress produces the drag force per unit area on a 
body traveling through the fluid. The shear stress at the surface provides a natural 
scale for velocities and distances near the wall for both laminar and turbulent flow. 
The friction velocity is deflned as 



/To 
P 



and the length scale, one "wall unit," is defined as 



= Iwu = — (2.4) 



The boundary layer thickness grows as the fiuid progresses over the body. For 
flow over a flat plate, the laminar boundary layer displacement thickness grows as 

5. ^ 1-72^^ 

The corresponding boundary layer Reynolds number, Re^ = increases with the 
boundary layer thickness until it reaches a critical value of Re^ ~ 420, at which 
point the laminar flow becomes unstable. This occurs at the corresponding critical 
boundary layer thickness of 

5^ ~ 420^ 
U 

and, for a flat plate, at the position 

when the boundary layer becomes unstable to the growth of TS waves. The Reynolds 
number describing the transition point is also often give as Re^ — corresponding 
to the distance from the plate edge at which transition occurs. 

Because the shear stress, tq, decreases as the boundary layer thickness, 6gg, grows, 
the wall units of the laminar boundary layer do not grow as quickly in x as the 
boundary layer thickness. The wall units for the laminar proflle grow as 

At the transition point these wall units are yl ~ 27^. 
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2.2 The boundary layer TS instability 



The linear stability analysis of laminar boundary layer flow over a flat plate was first 
calculated by Tollmien in 1929, Schlichting in 1933, and confirmed experimentally by 
Schubauer and Skramstadt in 1947. The stability analysis proceeds by considering 
an infinitesimal, z independent perturbation, tZt, to the laminar velocity field, wj^, of 
the form 

and solving the resulting eigenvalue problem for the wavenumber, frequency, a;, 
and perturbing profile, Ui{y). For a laminar velocity profile, Uq, determined from a 
solution of the Navier-Stokes equations, the imaginary component of oo first becomes 
positive when the boundary layer has grown to the thickness 5'^ ~ 4:20 jj, implying 
the onset of instability. The frequency and wavelength of the unstable TS waves at 
this onset of instability are found to be ^(to^) ~ -15^ and = || ~ 175^- The TS 
waves grow very rapidly beyond the transition point, and the boundary layer takes on 
a qualitatively different form. The new flow, now with large scale periodic variations 
in X, becomes unstable to a three-dimensional instability with periodic variations in 
z. The three-dimensional instability, best understood through a study of vorticity 
dynamics, evolves directly to turbulent flow. 

2.3 The evolution of boundary layer vorticity 

The evolution of the boundary layer may be more deeply understood via a description 
of the dynamics of vorticity. The vorticity held of the fluid is defined as the curl of 
the velocity. 

Taking the curl of the Navier-Stokes equation, (|2.1j ), and using a vector identity 
produces the vorticity equation for the incompressible fiuid, 

dtTJ- = dtZt + (l? ■ lo^ = ( • ^) 1? + uV^uf (2.5) 

The laminar fiuid fiow over the fiat plate is described simply as the creation of con- 
centrated vorticity at the leading edge and its diffusion, advection, and stretching 
downstream. 

2.3.1 Formation of vortex filaments 

The advected vorticity initially coats the fiat plate with a smooth layer of vorticity. 
The entire laminar boundary layer may be considered a vortex sheet of constant 
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strength, with sheet thickness increasing with x. From this point of view, the TS 
instability of the boundary layer is seen as similar to the Kelvin-Helmholtz instability 
of a shear layer. 

With the onset of the TS instability, the smooth vorticity of the sheet condenses 
into areas of higher vorticity separated periodically in x. These concentrations of 
vorticity then continue via nonlinear evolution to wrap and condense into spanwise 
oriented vortex filaments, evolution similar to the nonlinear development of a shear 
layer ||^]. This process of nonlinear development into vortex filaments is highly infiu- 
enced by small irregularities in the initial fiuid velocity field which produce significant 
variations from perfectly uniform vortex filaments. Although the shape and distri- 
bution may deviate significantly from the ideal, these filaments and their evolution 
represent the dominant structure of the boundary layer after the transition point, and 
play the central role in further boundary layer development. 

2.3.2 Three-dimensional vortex filament instability 

The initial vortex filaments, evolving from the vortex sheet after the TS instability, 
are separated by the distance A^. ~ 175^ in x, and are located at the initial height yo 
above the surface. A calculation of the three dimensional instability will show that 
this height must be equal to yo ~ 6wu. Since 3> yo, the effect of initial neighboring 
filaments on evolution is small and may be ignored |I0| . 

A straight vortex filament in the presence of a boundary is unstable due to the 
velocity field induced by its interaction with the surface. The filament velocity field 
is calculated by considering it to arise from an image vortex beneath the fiat surface. 
A plot of the image induced velocity field surrounding the filament, relative to the 
motion of the filament, (Figure shows the nature and direction of the stretching 
fiow. 

The filament is linearly unstable to deformations in this stretching plane, oriented 
at an angle of approximately 45° to the boundary surface. This three dimensional 
vortex filament instability is treated in detail in Chapter |^, and the effects of a wavy 
boundary are considered in Chapter ||. The filament over a fiat boundary is found to 
be maximally unstable to deformations of wavelength A^ — ITt/o, and is sketched in 
Figure |]|. 

A numerical investigation confirms that a filament with a small initial perturbation 
rapidly deforms into a sinusoid of wavelength A^ inclined at an angle of approximately 
43°. The expanding sinusoidal filament then leaves the realm of linear evolution as 
the bottom of the filament nears the surface and adheres to it, forming streamwise 
oriented vortex "legs," while the top of the filament continues to rise with linear rather 
then exponential growth as it advects downstream. In this manner a "horseshoe" or 
"hairpin" vortex is formed. Although the existence of hairpin vortices was postulated 
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Figure 2.2: The velocity field induced by the image vortex, relative to the motion 
of a straight vortex filament located at x = 0, y = 1. The image vortex is located 
at X = 0, y — —1. The self-induced velocity field and background fiow field are not 
included. 
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Figure 2.3: A plot of the maximally unstable deformation mode of a vortex filament 
over a flat surface. 



by Theodorsen in 1952 |T^, and has been examined in great detail in experiments by 
Head and Bandyopadhyay [l^ and others jl^, the vortex filament instability mecha- 



nism, the Crow instability ||Tl|, leading to the formation of these hairpin vortices, has 
only now been established ||10[- The linear Crow instability of a vortex filament pair 
was initially calculated, and confirmed numerically ||l^], for the case of two aircraft 
trailing vortices. The same analysis applies for the near wall filament, and is described 
in Chapter |^, with the second of the vortex pair replaced by the image filament. 
Recent experimental investigations indicate that autogeneration occurs in the fur- 



ther evolution of these vortex filaments ||T^. In this process, the legs of the hairpin 



vortex pinch together and form another filament, which then grows via the same in- 
stability, to align itself slightly behind the original. This process repeats and produces 
a long chain of hairpin vortices in line behind each original hairpin. These trains of 
tightly packed hairpin vortices, "hairpin packets," are the dominant structures in the 
turbulent boundary layer. 



2.4 The turbulent boundary layer 

The fully turbulent boundary layer over a fiat plate is a tumultuous collection of 
hairpin vortex packets and filament remnants. Since boundary layer evolution is 
highly nonlinear and sensitively dependent on initial conditions, any simple model, 
such as the one described with the idealized vortex filament picture, can only represent 
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the dominant processes of the dynamics and cannot capture the fine structure of small 
scale motions. Rather, a statistical analysis of fluid flow is employed to understand 
the processes of energy transport and induced surface drag. The effects of coherent 
structures, such as collections of hairpin vortices, may then be considered on the basis 
of how they effect the mean properties of the flow. 

2.4.1 Reynolds decomposition 

The mean velocity fleld in the turbulent boundary layer is found by taking the infinite 
time average of the velocity field 

1 /-to 

u = lim — / u dt 

io^oo to Jo 

The Reynolds decomposition is performed by decomposing the variables into the time 
average quantity and the remaining fluctuating quantity, u' , 

it = Tt +il 

This decomposition allows the steady, time independent, mean flow to be considered 
separately from the more complex underlying turbulent flow. The mean equations of 
motion are obtained by taking the time average of the incompressibility condition, 
( ^^ , and Navier-Stokes equation, (|2.1J ), to get, in component notation, 

d.Wi = (2.6) 
dtUi = UjdjUi = — dip + udjdjUi — dj (w^^ j (2.7) 



The last term, —dj [u'iUjj , the Reynolds stress term, represents the effect of the 
turbulent eddies on the mean flow. 

The mean flow of the boundary layer develops a proflle similar to the laminar 
flow proflle. A standard model of the mean velocity proflle for turbulent boundary 
layer flow has emerged from many decades of experimental research. The profile, 
function of y, breaks into several regions. The "inner region" extends to 
approximately 5wu, and is fitted by the curve ih;{y < 5wu) ~ m^,^, with ^ equal 
to the non-dimensionalized distance from the surface in wall units as defined by 
(^]4|). The "log layer" begins at approximately 30wu, and is fitted by the curve 

~ tt^ ^J^i-^) + 5). These curves are matched to each other in the "buffer layer" 
between 5 and 30wu. The log layer extends out into the "outer region," where it 
asymptotically approaches the free stream flow, U, (Figure [2^). The mean boundary 
layer thickness grows linearly with x as ^99 ^ •4^x, faster then the square root 
dependence of the laminar boundary layer. 
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Figure 2.4: The turbulent boundary layer mean velocity, u^{y), showing the various 
regions of the boundary layer. The mean velocity, along the horizontal axis, is shown 
as a function of the height above the surface in wall units, y+ = ^, along the vertical 
axis. 
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Unlike the laminar flow case, in which the primary energy loss is due to viscosity, 
the mean flow of the turbulent boundary layer is influenced more strongly by the 
energy loss to the underlying turbulent flow due to the Reynolds stress. In regions 
far from the boundary the effects of viscosity on the mean equations of motion are 
negligible, and the viscosity term may be ignored. 

2.4.2 Energy flow 

The mean energy transport equation is obtained by multiplying ([2?^ ) by Ui to get 



dt [ -UiUi] = di I — pui + 2vUjEij — u'{u'jUj — 2uEijEij + u'^u'^djUi 



in which the mean strain rate is Eij = \ {diU] + djut). This equation is best under- 
stood by considering it a description of the mean energy balance of an inflnitesimal 
fluid element. The change in kinetic energy of the element, dt QuTul), is equal to 
the energy transported across its boundary (the divergence term), plus the loss to 
viscosity, —2uEijEij, plus the loss to turbulence, u'^u'^djUl. 

The most significant mean strain in the boundary layer is dyU^. Hence, the dom- 
inant energy transport from mean to turbulent motion is due to the term u'^u'ydyll^, 
and large energy losses to turbulence are associated with u'^Uy < events in a region 
of large mean shear. This may be understood physically as the creation of turbulent 
flow by the lifting of low speed fluid near the boundary layer up and against the 
higher speed fluid in the upper boundary layer, a "low speed streak" and "turbulent 
burst," and by the converse, a "sweep," in which high speed fluid descends into the 
lower boundary layer ||14|| . This intermittent turbulent bursting and sweeping phe- 
nomenon is experimentally observed in the turbulent boundary layer, and measured 
to account for the majority, 60 — 80% of the turbulent energy transfer. These 

phenomenon are consistent with the hairpin vortex model of turbulent structures, 
and the large energy transfer may be considered as flowing from the mean shear into 
the stretching of the inclined hairpin vortex filaments in the boundary layer. 

The energy imparted to the turbulent development of the fiuid, that is, the growth 
of hairpin packets, must ultimately arise from the surface via surface drag. The mean 
surface stress is = /j,^ The mean velocity profile, u^{x,y), of turbulent 

boundary layer fiow is strongly affected by the Reynolds stress. The addition of 
the Reynolds stress term with u'^u'y < in the mean momentum equation, (|2.7|) , 
generates a much steeper profile than for laminar fiow, dyll^ ^ (^y^xiaminart ^^'^ hence 
much greater surface drag. 
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0.02 




Figure 2.5: The Reynolds stress, function of x along the line y = 1, 

z = threading a sinusoidal filament of unit circulation centered at y = 1, x = 
and oriented at 45°. The example filament has amplitude Xa = .75, ya = .75 and a 
wavenumber of k = .5. The induced velocity of the image filament is also included. 

2.4.3 Hairpin vortex contribution to Reynolds stress 

It is important to note that the naturally occurring hairpin vortices over a flat sur- 
face, inclined at an approximately 45° angle, produce the largest contribution to the 
Reynolds stress. The legs of such a vortex, oriented in the ±{x + y) direction, have 
large corresponding fluid velocities in the ±{x — y) direction, producing u'^Uy < 0. 
A plot of u'^Uy along a line threading an example vortex, (Figure [275|) , shows this 
explicitly. The legs of a hairpin vortex, inclined at an an angle 9, will produce a 
Reynolds stress, 

<^ ~ (- sin(0)) (cos(e)) = sin(2^?) (2.8) 

proportional to sin(26'). 

The velocity held around an advecting hairpin vortex, or vortex packet, directly 
explains the observation of streaks, bursts, and sweeps. The streak and burst events 
are described by the upwards and backwards flow beneath the vortex fllament heads, 
while the sweep events are associated with the passing of the vortex legs near the 
boundary. Furthermore, the well documented spanwise spacing of the streaks, A ~ 
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lOOwu ^7j, is directly explained as the wavelength of the maximally unstable mode 
of the vortex filament instability that arises from a filament born at yo ~ 6wu, the 
top of the high shear inner region of the boundary layer profile. 

A coherent model of turbulent motions over a fiat surface is formed by the con- 
sideration of the birth and development of hairpin vortices. The spanwise vortices 
condense during the TS instability, and rapidly develop into hairpins, via the three 
dimensional instability, to constitute a turbulent boundary layer. Further spanwise 
vortices condense in the high shear region, yo — Qwu, via autogeneration, and evolve, 
once again via the three dimensional instability, to form hairpin packets and fill the 
evolving boundary layer. 
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Chapter 3 

Filament evolution over a flat 
boundary 

The dominant structures in the turbulent boundary layer, the hairpin vortices and 
hairpin packets, develop via the three dimensional instability of concentrated patches 
of spanwise oriented vorticity. The vortex filament model of the developing turbulent 
boundary layer provides an elegant framework in which to calculate the specifics of 
this instability. 

The mean ambient vorticity at the high shear, inner region of the turbulent bound- 
ary layer, 

CO ~ —ZOyUx ^ —z — 

' y* 

accumulates into the thin, spanwise filaments centered at the filament height, yo, 
above the wall. The vorticity of each filament is considered to condense from a 
square patch of side yo centered around the filament center, producing a vortex with 
circulation 

r ^ — 2/0 
y* 

which propagates in the surrounding, locally constant, mean background fiow, u'' ~ 
u^iyo) ~ ill the x direction. 

In the region of the turbulent boundary layer dominated by vortical structures, 
above several wall units, the effects of viscosity are small compared to the dynamical 
interaction of the vortices with the wall and with themselves. The simplifying assump- 
tion is thus made that the filament propagates in an inviscid background, = 0, with 
the no slip condition at the wall relaxed to the free slip condition, Uy{y = 0) = 0, or, 
more generally, 

n • = (3.1) 

where n is the normal at the boundary, B. 
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3.1 The vortex filament 



Consider a single, nearly straight, infinite vortex filament, parallel to the z-axis, 
located in space above the flat surface at the position 

l{z,t) = {X{z,t),Y{z,t),z) (3.2) 

This time, t, and spanwise coordinate, z, dependent description of the filament posi- 
tion is adequate for the linear stability analysis in which the filament is not allowed 
to loop back on itself in the z direction. 

Ideally, the filament is represented as an infinitely concentrated region of vorticity, 

u^ = r dz'6\^ - R{z',t))^ (3.3) 

J ~oo CiZ 

The actual vorticity distribution is not singular, but is characterized by a finite vortex 
core size. This must be considered in calculating the self-interaction of the filament. 
However, the above approximation is valid in the case of small vortex core sizes, and 
the results of the stability calculation will display only a weak dependence on core 
size. 



3.2 Equations of motion 

The motion of the filament is obtained by using (|0|) in the vorticity equation, (|2.5|) , 
and dropping the viscous term to get 

dtl = dtl + u,\^ d,'^ = (3.4) 

— t — t 

in which the velocity field, u = w + u , which includes the background fiow and the 
filament induced fiow compatible with the vorticity distribution, is being evaluated 
at the vortex filament position, R{z,t). Inserting (p.2|) into (|3.4|) gives the equations 
of motion for the X{z, t) and Y{z, t) coordinates of the filament, 

dtX = -u^ d^X + 

dtY = -u,d,Y + Uy ^''•^^ 

The fiow induced by the vortex filament, 1?, must be calculated at each moment from 
the geometry of the filament and boundary. 
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3.3 Induced flow 



The filament compatible velocity field, , must satisfy u-^ = ^ x u-^ , with uj-^ given 
by (^.31), as well as satisfying the free-slip boundary conditions, (|3.1|). Since the 
fiuid is considered incompressible, the velocity may be written as the curl of a stream 

function, 1? = V x \1/ , with the stream function , \1/ = "^-^ + '^, consisting of terms for 
the filament induced fiow and the background flow, and background streamfunction 

= vf'yz for the current case. This gives, via a vector identity, the relationship 
between the vorticity and stream function. 



^ X X ^ = -V^^ + ^ (^"^ ■ ^ 



The gauge choice may be made without affecting the flow, . Finding the 

stream function, and hence the induced velocity, is then reduced to the problem of 

solving the vector Poisson's equation, V^'^^ = —UJ^-, with the appropriate boundary 
conditions. This solution may be written as 



by using the (dyadic) Green function satisfying 



-5'-^ {It - x')5i 



and the boundary condition 



V{z,/} 



(3.6) 



(3.7) 



(3.^ 



equivalent to it x G 



B 



0. 



Green's function satisfying ( |3.7| ) and ( |3.8| ) for a flat boundary at ?/ = is 

G^{{x,y,z),{x\y\z')) = 



(x-x'f+{y+y'f+{z-z'f 



1 

1 

1 

1 
0-10 
1 



(3.9) 



The construction of this Green function may be understood via the method of images. 
For each vortical element that acts as a vector potential source through the first term 
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in ( |3.9| ), there is a mirror image source beneath the surface that acts through the 
second term in ( p.9|) . The minus sign in Gyy arises because the mirrored element 
points in the opposite direction in y. The image ensures that the boundary conditions, 

( p.8|) , are satisfied, since it cancels out the components of parallel to the surface. 

The velocity field induced by the vortex filament and its image, as well as the 
background fiow, is found by taking the curl of the stream function to get 



u 



47r J I ^_ ^v^/ 



dR 
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X 



X 



Wiz') 



dW 
~dz' 



(3.10) 



which is recognizeable as the Bio-Savart law for the vorticity/current source, with the 
image vortex position defined as R^{z) = {X{z), —Y{z), z). 



3.4 Two dimensional motion 

The above equations simplify dramatically for the case of a straight filament. In this 
case, the filament position is independent of z, and can be written as 

%{t) = (xo (t),2/o (t)) 



The equations of motion, (|3.5| ), become simply 



dtXQ = 

dtVo = Uy 



and the induced velocity may be found by integrating (|3.10|) or solving Poisson's 
equation in two dimensions, and differentiating, to get 

_ „,f) r {y-yo) , r {y+yo) 



(x-xq) r (x-xq) 

27r (^x-xoY+(y-yoY 27r [x-xoY +(y+yoY 



Evaluating this induced velocity at the filament in the equations of motion gives 

dtVQ = 



which implies that a straight filament travels over the fiat surface at a constant height, 

1 



with a constant velocity of m** + in the x direction. 
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3.5 Hamiltonian formulation 



It is interesting and very useful that the dynamics of the vortex filament motion, 
(^■5|), have a Hamiltonian formulation ||2T]|. Taking the dynamical variables to be 
X{z,t) and Y{z,t) and the canonical, instantaneous, fundamental Poisson bracket to 
be 

{Xiz),Yiz')} = biz 



the equations of motion are obtained as 

{X,H} - 



dtX 
dtY 



V-^ = -u,d,X + u, 



r SY 

1 6H 

'v &x 



(3.12) 



(3.13) 



in which the Hamiltonian is the kinetic energy of the filament induced flow as a 
functional of X and Y , 



H 



y d^x xP\2 (^^ xpy X 

U d^x (p + 2^^ ■P'^ 

11 J J dzdz'^G,, (tiz) , -niz')) + r / dz^ ■ ^ 



(3.14) 



5H 



Note the common abuse of notation by which a variational derivative such as gx{z) 
deflned implicitly such that 



IS 



SH 



6H 



The existence of Hamiltonian dynamics for the motion of the filament will facilitate 
the linear stability analysis as well as aid in the calculation of the equations of motion 
in new variables via canonical transformations. 



3.6 Self interaction 

The equations for the three dimensional motion of the filament, (|3.4|) , are presently ill 
defined because the induced velocity, (|3.10|) , diverges at the filament. This is remedied 
by using a filament with a finite core size. Rather then make a smoothing correction 
to the vorticity distribution, ( WM , with the associated complication of dealing with 
core dynamics ||22|||2^, a cutoff is imposed on the self-induction term for the stream 
function § 



dRj 
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where the [S] stands for the region z' < z — aS and z' > z + aS, a is the vortex core 
radius, and (5 is a numerical factor. The approximate value of 6 is determined by 
establishing consistency with the known self induced velocity of a solid cored vortex 
ring, and is found to be 5 ~ ^e*. 

For small core sizes this cutoff method is equivalent to the modification of the 
Green function to H 



Airy {x~x'Y+{y-y'Y + {z-z'Y+e~^a 



^{x-x'f+{y+y'f+(z-z'? 



1 
1 
1 
10 0' 
0-10 
1 



(3.15) 



The use of this approximate, non-singular Green function is superior to the cutoff 
method for calculations. It is used in (|3.14|) to form the new, non-singular Hamilto- 
nian. 

The radius of the vortex cores, a, for the initial spanwise filaments must be es- 
tablished by direct experimental observation. An inspection of slides of turbulent 
boundary layer fiow indicates core sizes on the order of a ~ Iwu [l^. For a filament 



located at an initial height of ~ Qwu this corresponds to a relative core radius of 



a 

yo 



1 
6 



(3.16) 



Although this quantity is only roughly estimated, the filament instability of inter- 
est, the instability induced by the filament interaction with the wall, is only weakly 
dependent on core size. 



3.7 Linearization 

The equations of motion may be linearized for small perturbations from the straight 
filament solution. An analysis of the filament behavior for small perturbations will 
display the nature of the three dimensional instability, including its angle and growth 
as a function of spanwise wavelength. This linear analysis will later be confirmed and 



extended by a fully nonlinear numerical simulation in Section |3.10 . 

The equations of motion, (|3.5| ), are linearized around the two dimensional filament 
motion to get, 

X{z,t) = xo{t) + ex{z,t) 

Yiz,t) = yoit)+Eyiz,t) ^■'■''> 
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where e is a bookkeeping variable used to keep track of the order in the small per- 
turbative variables, x and y. The analysis of linearized motion is facilitated by the 
Hamiltonian formulation. The Hamiltonian, (|3.14 ), for the filament over the flat 



boundary is expanded order by order in e to obtain 

in which if^*^^^ is the part of the Hamiltonian quadratic in x and and H^'^^^ is zero. 
The equations of motion, ( p.l3| ), are then similarly expanded order by order in e. The 
zeroth order equations are the two dimensional equations of motion, ( 3.11j ). The first 



order equations produce the linearized equations of motion for the filament. 



d,x = i^^"'^' 
= -r Si 



1%F(2) (3.19) 



which depend only on the quadratic part of the Hamiltonian, 

H^i^) = -[ [ dzdz' ( + G^^A (3.20) 



2 7 7 V 9^ "''^ 9z' 

in which ^ = {x, y), Greek indices {a, (3} range only over {x, y}, and G^j^^ and Gfj^^ 
are the modified Green function, ( p.l5| ), components of zeroth and quadratic order in 

e. 



Glf^ = G^^0{z),^{z')) 
GT = li.GU'^iz),t{z')) 



£=0 



e=0 



3.8 Diagonalization 

The linearized equations of motion, ( |3.19|) , along with the quadratic Hamiltonian, 
( p.20| ), comprise a set of linear, integro-partial differential equations for the dynamical 



variables x{z, t) and y{z, t). This set of equations is reduced to a set of linear, ordinary 
differential equations by applying a Fourier transform in z to obtain the new complex 
dynamical variables, x{k,t), defined as 

1 f°° 

x{k, ^) = y d,z e~'^''x{z, t) 

with 



oo 

ikz' 



x{z,t) = / dke'^'x{k,t) 



oo 
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and y{k,t) defined similarly. This is a canonical transformation. The new fundamen- 
tal Poisson bracket becomes 



1 r(f^/^£W£CHl - Sy*{k')Sx{k) \ 
r J !_ Sx{z) Sy{z) Sx{z) Sy{z) J 



(3.21) 



The next step is to show how the quadratic Hamiltonian, ( p. 20]) , is diagonalized 
by this transformation, and calculate the resulting equations of motion. Since the 
Green function, ( p.l5|) , within the Hamiltonian is explicitly dependent only on z — z' 
and not on z + z', it is useful to write it as a function oi ^ = z — z' and perform the 
Fourier transform in this variable to obtain 

G^ix,y,C,x',y') = [ dke'^^G^{x,yXx\y') 



in which the Fourier transformed Green function is calculated as 

G^{x,y,k,x\y') = ^ j die-^^^G^ {x,y,i,x\y') 



^Ko{k^{x - x'Y + {y- y'Y + e'^a^ 



(2^) 



^K,{k^{x-x'f + {y + y'f) 



1 
1 
1 
1 
0-10 
1 



with Kq the modified Bessel function. The quadratic Hamiltonian is then written as 

^Fi2) ^ J J J dzdz'dke^'^^ S^^G^J\x{z),Y{z),k,X{z'),Y{z')f-^ 
+G^P{X{z),Y{z),k,X{z'),Y{z'))} 

The and G^^^^ are computed, the Fourier integrals are substituted for x{z) and 

y{z), and the resulting quintiple integral vanishes in a poof of delta functions to give 
the diagonalized Hamiltonian in the Fourier variables, 

p2 „ 

i7^(2) = _ J dkk^ {A{k) x*ik, t) X (k, t) + B{k) y*{k, t) y{k, t)} (3.22) 



with 



^ = S{ka) + {-^ - Ko{2kyo) - ^K,{2kyo)} 



B 



S{ka) + {- 



2i^^i(2%o)} 



(3.23) 



and the vortex filament self interaction term 



S{ka) 



- , } + Ko{e—ika) + -^Ki{e-ika) 

(e-^ka) e ^ka 

0.183 - I ln{ka) + 0{{kaf) 
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(3.24) 



The linear equations of motion for each mode are then easily computed from ( p.22|) 
and (|3^ ) to be 



dtx{k,t) 
dty{k, t) 



1 srpJ^) 

2nr Sy* 

1_SH^ 

2nr Sx* 



(3.25) 



3.9 Linear instability 



The set of linear ordinary differential equations of motion for each mode, ( p. 25]) , are 
reduced to an eigenvalue problem by assuming a solution of the form 



X 

y 



xie' 



crt 



crt 



to get 



a 



Xi 









The two resulting eigenvalues are 



± 



Tk^ 



27T 







V-AB 



Xi 

m 



(3.26) 



with the angle, up from x, of the corresponding eigenvectors equal to 



arctan(-^) 

Xi 



arctanf 



r y/^AB . 

iri B ' 



All 
with 



The positive eigenvalue, or growth parameter, (T+, is plotted in Figure ^TT 
small wavenumber perturbations of the filament are unstable up to ky^ ~ .53, 
the maximum instability occurring at ky^ ~ .36. For wavenumbers above .53 the 
eigenvalues, (p.26|) , become imaginary, corresponding to filament perturbations that 
rotate rather then grow. The unstable eigenmodes corresponding to (T+ grow at an 
angle, 6*+, plotted in Figure All unstable modes grow in planes with angles of 0° 
to 90° up from x, with the maximally unstable mode growing at an angle of 6^+ ~ 43°. 
The wavenumber of the maximally unstable mode is dependent on the choice of core 
size, ( p.l6|) , for the calculation. However, the instability is only weakly dependent on 
this choice. For example, a core of only one tenth the current estimated size of ^ = | 
produces a similar instability plot with the maximum at ky^ ~ 0.28. Also, the angle 
of the maximally unstable mode remains 6*+ ~ 43° over several orders of magnitude 
of core size. 
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Figure 3.1: The growth parameter, (T+, in units of p^, as a function of spanwise 
wavenumber, kt/o. The maximum at kyo ~ 0.36 is marked. The core size was set to 

be — = I to obtain this plot. 
2/0 6 ^ 
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e 




kyo 

Figure 3.2: The instability angle, 6^, as a function of spanwise wavenumber, ki/Q. 

The maximally unstable wavenumber, kyo ~ 0.36, is marked, indicating the angle of 

the maximally unstable mode to be 9+ ~ 43°. The core size was set to be — = 4 to 

^ yo o 

obtain this plot. 
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For any small initial perturbation of the straight filament, the high wavenumber 
modes rotate about the core while the low wavenumber modes expand in the unstable 
plane, 0^, and contract in the stable plane, 6'_. The maximally unstable mode, ki/Q ~ 
0.36, expands exponentially faster then the others, and rapidly becomes the dominant 
deformation of the filament. The filament evolves into a sinusoidal shape inclined at 



approximately 43°, with a wavelength of A = ~ 17|/o, as plotted in Figure ^ 



Although the linear analysis gives a precise description of vortex filament motion for 
small perturbations, the full equations of motion, ( |3.5| ), are highly nonlinear; and it 
is expected that nonlinear mechanisms will determine the dynamics after the initial 
perturbations grow to significant size. 



3.10 Numerical simulation and nonlinear evolution 



Numerical simulation provides an effective means of exploring nonlinear behavior, 
as well as confirming linear analysis. The fully nonlinear equations of motion for a 
vortex filament lend themselves readily to numerical analysis. However, the equation 



for vortex position, (|3.2|) , must be modified to a fully parametric representation to 
allow for the possibility of the filament doubling back in the z direction, 

t{s,t) = {X{s,t),Y{s,t),Z{s,t)) 

in which s is an arbitrary parameter along the filament. The equation of motion for 
the filament is then simply 

dtR{s,t) = Ttl^ 



with the velocity field given by the Bio-Savart integral, (|3.10|) , with a finite core size. 
The full set of nonlinear integro-partial differential equations is 



dtR{s,t) 



R{s,t)- R{s',t) 

'M{s,t)-R^{s',t) 



3 „ 



X ds'R{s',t) 



R{s,t)-R {s',t) 



X ds'R^{s',t) 



with the image position R {s,t) = {X{s,t), —Y{s,t), Z{s,t)). These equations are 
then discretized in s and integrated forward in time. Care is taken to insert new 
discretization points along the filament as it is stretched, to preserve the accuracy of 
the scheme. 
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Figure 3.3: The vortex filament evolving nonlinearly from a small initial Gaussian 



perturbation. This figure may be compared directly with Figure 2.3. 



An example of the resulting filament motion is shown in Figure p.3| . For this case, 
the initial filament was taken to have a small Gaussian perturbation in the x direc- 
tion, i?(s, 0) = ^.le"^*) , to simulate the introduction of a small disturbance 

to the fiow. This Gaussian has a small projection onto the maximally unstable mode 
of wavelength A ^ 17?/o, which rapidly grows to dominate the filament geometry, con- 
firming the accuracy of the linear analysis. The nonlinear effects begin to contribute 
significantly to the evolution as the legs of the filament approach the surface. The 
legs are drawn out in the —x direction, producing counter rotating streamwise vor- 
tices near the wall. And the initially planar sinusoidal curve of the filament becomes 
slightly bowed as it expands. 

Although this numerical simulation gives a satisfactory description of filament mo- 
tion for short times, it is insufficient to address long time evolution questions such as 
the details of filament autogeneration and multiple filament interaction. A simulation 
of autogeneration along these lines would require more detailed core dynamics as well 
as the use of vortex filament surgery. Alternatively, a three dimensional simulation 
of the Eulerian vorticity field of an autogenerating hairpin produces excellent results 
p!3|| . However, the present numerical model serves as a successful confirmation of the 
linear result, displaying the dominant evolution of 6*+ ~ 43" modes with wavelength 
A ~ ITyo, and provides a consistent picture of the evolution of single hairpin vortex 
features. 
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3.11 Comparison with coherent structures in the tur- 
bulent boundary layer 

The qualitative agreement between the analytic results and the reported experimen- 
tal observations of hairpin vortices leaves little doubt that this filament instability is 
responsible for the development of hairpin vortices in the turbulent boundary layer. 
The observed 45° inclination is given a solid mathematical foundation, and the nu- 
merical simulation of the evolving filament, |3]^, produces a geometry that conforms 
excellently with the reports of observers and conceptual models as reviewed by Robin- 
son ||19|, and more recently elucidated by Zhong and others p3, and by Delo and 



Smits ||5|. 



The analytic calculation of the existence of the maximally unstable wavenumber, 
kyo ~ .36, together with the experimental observation of spanwise wavelengths of 
lOOwu for these structures in the boundary layer, demands that the average initial 
height of the vortex filaments be equal to yo ~ 6wu. This is reasonably located in a 
region of high mean shear in the turbulent boundary layer, and agrees well with the 
measured location of maximal turbulent energy production |26|| . 

However, the true power of the newly established mathematical model for vortex 
filament dynamics in the boundary layer lies in its potential for describing and pre- 
dicting new phenomenon. In the next chapter the existant mathematical framework 
will be extended to the case of a corrugated boundary surface, and new predictions 
will be extracted regarding the subsequent dynamical behavior of the vortex filament. 
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Chapter 4 

Filament evolution over a wavy 
boundary 



Several decades of intense academic study have provided an extensive 
knowledge base of the simple canonical case. The immediate need is to 
learn to utilize this store of information in the context of boundary-layer 
modeling and control methodologies, with the eventual goal of practical 
application to engineering problems involving real-world, noncanonical 



boundary layers. - Stephen K. Robinson |19 



In the last paragraph of his excellent review article on coherent motions in the turbu- 
lent boundary layer, Stephen K. Robinson urged researchers to build upon the foun- 
dation of knowledge in vortex filament behavior and to forge ahead into the uncharted 
territory of control methodologies. This chapter addresses one such methodology: the 
use of a static, wavy boundary to control vortex filament evolution. 

The study of this control methodology is motivated by the existence of such sur- 
face corrugations on the bodies of dolphins, with the physical parameters of back- 
ground fiow speed and surface ridge geometry taken directly from measurements of 
these animals. A study of vortex filament dynamics above such a boundary, using 
the tools developed in Chapter H, elucidates the effects on turbulent boundary layer 
development and the potential for drag reduction. 

A vortex filament propagating over a wavy surface is subjected to rapid periodic 
forcing due to its interaction with the surface. The filament is near the surface over 
the ridge peaks and far from the surface over the ridge troughs. Also, the unstable 



manifold, shown in Figure oscillates in direction and magnitude as the surface 
beneath the filament is inclined to the horizontal. This rapid parametric forcing is 
similar to the pondermotive forcing of a charged particle in an electromagnetic wave, 
and dynamically similar to the parametric forcing of an inverted pendulum p7[|,||28||. 
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The pondermotive forcing has a dramatic effect on the instability, as it does in the 
case of the pendulum. 



4.1 The boundary and the vortex filament 

The wavy boundary is taken to be an infinite surface, beneath the filament, 
uniform in the z direction, and periodic in the x direction. Hence, 

^(x, 2) = (x, B[x), z) 

with B{x) a smooth, periodic function with period L and trough to peak ridge height 
h. The surface normal is 

n = , - (—B'(x)x + y) 
and the free slip boundary condition, n ■ ^|-^ = 0, is assumed. 



The cutaneous ridges of Tursiops Truncatus, (Figure |1.2| ), are near sinusoidal, 
with slightly steeper ridges then troughs. The ridges of a live dolphin at rest are 
measured, 25cm behind the blowhole ||5[, to have an average height of /i ~ .015mm 
and an average spacing of L ~ .61mm, producing a height to length ratio of ^ ~ .025. 
It is expected that a dolphin will increase or decrease the ridge height, using cutaneous 



muscle, to accommodate changes in swimming speed, with the ridge height to length 

h 
L 



ratio possibly ranging up to 4 ~ 0.2. 



Dolphins swim at a cruising speed of 10^ ~ 3^, with a top speed of 50^ ~ 
14y. For a dolphin swimming at U = 3^, the friction velocity at, and slightly 
beyond, the transition point is ~.ll^. A wall unit at this point is thus 

?/* = ^ ~ .009mm, and the vortex filaments are assumed to form at a height of 
HQ ~ 6wu ~ .054mm. The vortices at this height propagate in a background fiow 
field of velocity 

U ~ ^ .00 — 

y* s 

and have a circulation of 




Tc^-—y^ = -u{ — ] ^-3.6x10"'"' 

y* 

Although these initial values for the filaments are calculated for a fiat plate, they will 
be taken as the approximate values for filaments evolving over a wavy boundary. 

The filament position as a function of z and t is, for the case of a wavy boundary, 
still written as 

'n{z,t) = {Xiz,t),Yiz,t),z) (4.1) 
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with the dynamical variables X{z,t) and Y{z,t) satisfying the equations of motion, 

(PD, 

dtX = -Uz dzX + 

dtY = -u,d,Y + Uy ^^■^> 

However, the induced velocity, as well as the background velocity field, differs from 
the field for a flat plate, (|3.10|) , because of the wavy boundary. 



4.2 Orthogonal curvilinear coordinates 

The existence of a wavy boundary significantly complicates the calculation of the 
induced velocity field over the surface. However, the problem is simplified by the 
introduction of orthogonal curvilinear coordinates that match the boundary. 

A coordinate transformation is made from x and y to new coordinates, u and f , 
such that 



x{u,v) = ^ arctan (tan tanh (t; + + 2Lp 

^ In (4 cos^ (^) cosh^ (f {v + vb)) (4.3) 
+4sin2(^) sinh^ (f (v + vb))) - vb 



y[u,v) = f 



in which the term with p, equal to the integer part of | (^^ + , is necessary to give 
a smoothly increasing x for u > L, and Vb is the ridge height parameter, related to 
the ridge height, h, by 

.^ = --ln(^tanh(-)j 

The coordinate lines of constant u and v, mapped into the xy plane, are shown 
in Figure ^T|. This coordinate transformation approximates the wavy boundary of 
dolphin skin, 

B{x{u,v))l^Q = y{u,v)l^Q 

with peaks slightly sharper then troughs. 

The transformation, (|4.3| ), is a conformal transformation derived from the complex 
analytic mapping, 

({w = u + iv) = X + iy = — In cos ^— {u + iv + ivs)^^ + 2Lp — ivs (4.4) 

and thus preserves the orthogonality of the coordinates. The great advantage of using 
a conformal transformation is the ease of calculating the two dimensional Laplace 
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Figure 4.1: The conformal coordinate mapping of lines of constant u and v into the 
xy plane, with a chosen height to length ratio of ^ = .2. 

operator, and solutions to the Laplace equation, in the new coordinates. Because the 
analytic function, (|4.4|) , satisfies the Cauchy-Riemann equations, 



duX = d^y 
dyX = -duV 

the two dimensional Laplace operator, in the new coordinates, becomes 

in which the function s^iu^v) is the scaling factor, or Jacobian, of the conformal 
transformation - a measure of the area change of the new curvilinear coordinate grid 
- and is equal to 

s^{u,v) = \C'iw)\^ = (dux) (dyy) - (duy) (d-ax) 

_ 2cosh(^(^+^B)) _ 

~ cosh ( ^ (»;+i7s ) ) +COS ( ) 



for the particular coordinate change, (|43|). 

Using the new coordinates, the background potential flow satisfying the free-slip 
boundary conditions over the wavy surface may be obtained from the simple vector 
potential, ^(m, v) = u^vz. This vector potential produces an irrotational background 
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velocity field, X ^f^, which satisfies the incompressibility condition, 

as well as the free-slip condition, n -v^ ^ = 0. The lines of constant v in the new 

coordinate system correspond to streamlines of this background fiow over the wavy 
surface. 

Although it is possible to continue to consider the motion of the filament in carte- 
sian coordinates above the wavy boundary, using the curvilinear coordinates only for 
the purpose of velocity field computation, the analysis is greatly simplified by trans- 
forming the dynamical variables of the filament, X and Y, into the new coordinate 
system as well. This obviates the need to convert back and forth between curvilinear 
and cartesian coordinates. 



4.3 Non-canonical transformation 

A non-canonical transformation is made from the old dynamical filament variables, 
X and Y, to the new variables, U{X,Y) and V{X,Y), via the inverse of the coordi- 
nate transformation, (|4.3|). The fundamental Poisson bracket in the new coordinates 
becomes 

iuyz),vyz)^ j j uz uz \sx(z") SY(z"') 5X{z") 5Y(z"') j r"y^ ^ ^ ^4 g^j 

with the scale factor, s^, appearing as a direct result of the compression and expansion 
of fiuid as it is carried over the wavy surface. The position of the vortex filament in 
uvz coordinates is now given by 

'niz,t) = iUiz,t),V{z,t),z)^^^ (4.6) 

and the Hamiltonian, (p.l4| ), in the new coordinates is written as 

with the indices now ranging over {u, v, 2}, and the new Green function for the 
wavy boundary. 

The boundary conditions for the wavy boundary Green function, 

^^,kn,GY^\^ = V{z,/} (4.8) 

take a particularly simple form in the curvilinear coordinates, in which n = v, with 
( [4.8| ) simplifying to 

41 



which is ^) X 
equation, 



0. The wavy Green function must also satisfy Poisson's 



D = 



V72 (^W 

^ ik^kj 



(4.10) 



with the three dimensional Laplace operator in curvilinear coordinates taking the 
form 



'^•) - X "^x 



1 
1 
1 



+ ^{{dus)d,-{d,s)du) 



1 
1 




(4.11) 



0-10 
1 




Although the solution to ( [1.10| ) is simple for the case of two dimensional motion, 
because the z dependence is integrated out, the full three dimensional solution is 
complicated by the position dependent Laplacian, and an approximation to this so- 
lution will be made in Section WM 



4.4 Two dimensional motion 

The goal of this analysis is to obtain a description of the behavior of a vortex filament 
over a wavy surface. The dynamics of a straight filament must be determined, and 
then extended via a small perturbation, to predict the filament behavior. The be- 
havior of the straight filament is determined by considering the motion in the proper 
curvilinear coordinates. A two dimensional Hamiltonian formulation, which arises 
as the limit of the full three dimensional Hamiltonian formulation, facilitates this 
analysis. 

The position of a straight filament, parallel to the corrugations of the surface, is 
described by 

^ (t) = («o (t) , ^0 (t)).. 

The dynamics of the filament are described via Hamiltonian dynamics, with the 
fundamental Poisson bracket, arising from (|4.5| ), equal to 
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and the Hamiltonian, from ([4.7|), equal to 

in which the two dimensional Green function for the wavy surface, 



is the solution to 



(4.12) 



'2 



y2(^iy = 1 1^ 92 d'^ \q\ 

= (4.13) 

and satisfies the boundary conditions, = 0. (Note the appearance of the scale 

1^0=0 



factor, ^, when the delta function in ( ^.13| ) is converted from a distribution over flat 



to curvilinear coordinates.) The resulting equations of motion for the straight filament 
propagating over the wavy boundary are 



dtvo = {vo,H^} = 



(4.14) 



The filament travels along a streamline of constant Vq, with Uo{t) steadily increas- 
ing with period T in time, UQ{nT+t) = nL + UQ{t). The first order, nonlinear ordinary 

differential equation for uo{t), ( |4.14|) , is integrated numerically to obtain the filament 

— ^ 

position as a function of time. The straight filament, R2 (t) = (mq (t) ,fo)„„, moving 

along the streamlines of constant Vq with varying velocity, (Figure |42|) , is used as the 

basis for small three dimensional perturbations about this solution. 

It is also possible to integrate the equations of motion, ( |4.14| ), analytically, rather 



then numerically, to obtain an equation in closed form for the filament postion, Uq, 
as a function of time. Carrying out the integration gives 



+ tanh(^(.o+..)) ^"'^^^ (^^^ (l"°J ^^^^ + 

which may be solved for Uo{t). However, solving this transcendental equation for Uo{t) 
is no less difficult then integrating the equations of motion numerically, and equation 
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Figure 4.2: The path of a straight vortex filament at height ^^o = 1 propagating over 
a wavy boundary of wavelength L = 10 and height to length ratio j = .2. 



( ^1.15| ) is hence used only to calculate the period of filament motion over the surface, 

L ( 2 



T 



-1 + 



+ 2^^) V tanh {v^ + VB) 
which is then used to set up the numerical integration of (|4.14| ) over one period. 



4.5 Linearized equations of three dimensional fila- 
ment motion 



The Hamiltonian, (^?7|), with the Poisson brackets, (|4.5|) , generates the dynamics of 
the vortex filament above the corrugated surface. 



dtU{z,t) 
dtV{z,t) 




1 



5H^ 



Ts'^{uy) sv 



5H^ 



(4.16) 



These equations of motion are linearized about the straight filament solution, as was 



done in Section |3.7| , with the breakup of U and V into two dimensional and small 
three dimensional terms, 



U{z^t) = uo{t) + eu{z^t) 
V{z,t) = vo + ev{z,t) 

The Hamiltonian is expanded order by order in e to obtain 



(4.17) 



(4.18) 



The equations of motion, (|4.16| ), are then similarly expanded order by order in e, 
with the first order equations producing the linearized equations of motion for the 
filament, 

Otu - -i^[^u+^v 

dtV = 



5u 



44 



in which and its derivatives are evaluated at the filament center, {uo,vo). The 
quadratic part of the Hamiltonian is equal to 



2 



dz 



dz' 



(4.20) 



in which Ft = {u,v), Greek indices {a, (3} range only over {u,v}, and G^f^^^ and 
GYJ"^^ are the modified Green function components of zeroth and quadratic order in 

GZ(ltiz),tiz')) 



e=0 



e=0 



4.6 The approximate Green function 

The calculation of the linear equations of motion for the vortex filament over the 
wavy boundary, (|4.19| ), requires the use of the Green function solution to the singular 
Poisson equation, ( |4.10| ), satisfying the boundary conditions, (^.91). This Green func- 
tion may not be easily found in closed form because of the non-trivial nature of the 



position dependent Laplace operator in the wavy coordinates, ( ^.11| ). However, it is 
possible to approximate this Green function, and to improve upon this approximation 
by employing a Neumann series. 

The approximate Green function is required to have the following properties: it 
must reduce to the fiat boundary Green function as ^ — >^ 0, it must integrate to 
the correct two dimensional Green function, (|4.12|) , and it must satisfy the boundary 
conditions at the surface, ( [1.91) . The fiat boundary Green function, (|3.15| ), with 
curvilinear coordinates as arguments. 



G 



w 



G^{{u,v,z),{u\v\z')) 



(u-u')'^+{v-v')'^ + {z-z')'^+e 



47n/ (u-u')'^ + {v+v'f + (z-z')'^ 



1 






1 





-1 







1 





1 



(4.21) 



satisfies these conditions. It is a good approximation to the exact Green function for 
the wavy boundary when applied in the case of small amplitude filament perturbations 



and small ^. 



Although the approximate Green function, (|4.21|) , is not an exact solution to the 
Poisson equation, it is an approximate solution, in that it is a solution to 



4v 



2F 



{{u, V, z), {u', v', z')) = --6iu - v!)b[v - v')b{z - /)4 



(4.22) 
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with the "flat" Laplace operator, V^^ = 8"^ + dl + dl. Equation (|4.22| ) approximates 
the Poisson equation, ([4.10|) , which may be written as 



= -j,Siu-v^6iv-v')Siz-z%, 
= —6^ {if — x')6ij 



(4.23) 



by breaking the Laplace operator in curvilinear coordinates, (|4.11|) , into a "flat" part, 
V^^, and "wavy" part. 



V 



2W 



S2y2 _ y2F 

is' - 1) dl 



1 
1 
1 



which goes to as ^ ^ 0. 
4.6.1 Neumann series 



1 
1 




(4.24) 
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If necessary, the approximate Green function, ( [4.21|) , may be improved by adding 
terms from a Neumann series expansion. By rewriting ( |4.23|) as 



W'^G^ = -6iu - u')6iv - v')Siz - z')6,, - = -A, («, v, z) 



(4.25) 



considering p to be the source for a "flat" Poisson equation, and using the "flat" 
Green function, , to "solve" for the G^ on the left hand side of (|4.25|) , a Fredholm 
equation of the second kind for is obtained. 



J J J du"dv"dz"Gf, {{u, V, z), {u", v", z")) p,,{u", v", z") 
Gf^+JJJdu"dv"dz"Gr,Vlf'G]^ 



This integral equation may be solved for G^ by recursive iteration, producing a 
Neumann series, 



G. 



w 



Gi,+ 



du"dv"dz"GiW^^Gf^ + .. 



(4.26) 
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which will converge for small V^^, and hence for small ^. 

In practice, the integrals in (|4.26|) cannot be written in closed form, and must 
be computed numerically to determine the corrections to the Green function. This 
procedure is not necessary to the calculation of the linear instability of the filament 
over a surface of small ^, which uses only the simple approximate Green function, 
~ iiu,v,z),{u',v',z')). 



4.7 Fourier transform and equations of motion 

The diagonalization of ([4 .201 ) is carried out via a Fourier transform in the z coordinate. 



as in Section with the small dynamical variables of the perturbation becoming 
the Fourier transformed pair, u{k,t) and v{k,t). The use of the simple approximate 
Green function produces an approximate quadratic Hamiltonian identical to (|3.22|) , 



HW{2) ^ ^ 1 (iA; P {A^{k) u*{k, t) u {k, t) + B^{k) v*{k, t) v{k, t)} 



with A and B given, as before, by 



with the self interaction term given by ( p.24| ) . 

The linear equations of motion, ( |4.19|) , for each mode then become, in matrix 
form, 



u 


1 


V 





a 

du 







2n 





u 




V 



(4.27) 



with and its derivatives evaluated at {uo{t),vo). The effect of the wavy boundary 



is apparent in the difference of (|4.27|) from the equations of motion in the presence 
of a flat boundary, (p.25|) . The most important difference is the appearance of time 



dependent terms in the matrix coefficients, due to the periodic motion of the straight 
vortex filament over the wavy surface, Uo{t). This affects the filament primarily 
through an overall time periodic amplitude stretching, ^2(^Q(t-) , due to the expansion 
and contraction of the flow streamlines over the surface. The perturbative mode 
amplitudes, u and v, are also affected independently by the changes in the flow 
velocity fleld surrounding the propagating fllament - a result of the velocity held, 
and hence the uv frame, rocking rapidly back and forth as the center of the filament 
travels over the undulations. 

Although the method of images is not accurate for a wavy surface, it is a good 
approximation when the filament height is small compared to the period, 4 <^ 1, 
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and provides an alternative view of the effects of the wavy surface on the filament 
instability. As the filament center travels over the surface, following a streamline, 



(Figure ^.21 ), it travels periodically closer to and farther from the boundary, with 
its image filament drawing correspondingly nearer and farther. This rapid, time 
periodic, change in the most fundamental parameter, the vortex filament height from 
the surface, has a significant effect on the evolution of the vortex filament perturbation 
amplitudes. 



4.8 Floquet analysis 



The linear equations of motion, ( [4.27| ), governing the evolution of the vortex filament 
perturbation for each wavenumber, k, contain coefficients that are periodic in time, 
with period T. Floquet's theorem implies that the general solution to these equations 
may be written as 



uit) 
v{k) 



ae 



ait 



.hit) 

f2{t) 



+ be 



U2t 



91 it) 

92 it) 



(4.28) 



with complex Floquet exponents, {<Ji,a2}, time periodic functions, {/i, /2, fi'i, fi'2}, 
and coefficients , {a,b}, determined by the initial conditions. 

An inspection of the time reversal symmetries of the coefficients in the equations 
of motion provides further information about the solution, allowing the second Flo- 
quet solution set, {(71, (72} and (72, to be related to the first. The symmetries of the 
coefficients in (|4.27] ) are 





-t) = 


s\t) 


dtuo 


-t ~ 








d 




du 


-t 


du 


d£_ 








dv 


-t 


dv 


t 



Hence, the existence of a solution with the time dependence 



u 




' fi (t) ' 


V 




J2{t) _ 



implies - through time reversal of the equations of motion - the existence of the 
solution 



u 

V 



'(Jit 



h i-t) 

-f2 i-t) 



Since the general solution. 



consists of only two independent terms, the second 
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term must accommodate this time mirrored solution, and hence 



gi{t) = M-t) 

92 (t) = -f2{-t) (4.29) 

It is not surprising that the Floquet exponents are inverses, since the motion is due 
to Hamiltonian flow and to coordinate frame rotation. 

The stability of the fllament is now completely described by the Floquet solution. 
Since the exponents are inverses, the perturbation has either neutral stability, with 
= 0, or a pair of unstable and stable modes, with exponents a± = ± |(Ti|and 
^{(Ji) = 0. A complete description of the linear behavior of the perturbation to the 
vortex fllament requires only the calculation of fi(t), /2(t), and cxi. 



4.9 Numerical implementation 

Although analytic methods exist to calculate the Floquet exponent, cxi, they are not 
practical for the current case because of the complicated time dependent coefflcients 
in the equations of motion, (|4.27|) . Rather, a numerical solution over one period, T, 
of the motion is used to determine the Floquet functions, /i and /2, and Floquet 
exponent, cti, as a function of wave number, k, fllament height, Vo, background flow, 
U, and surface geometry, h and L. 

Of primary interest are the size of the growth parameter, 0"+, if non-zero, and the 
angle, 9 = arctan(^), at which the perturbation grows over one period. To determine 
these quantities, the equations of motion, ( [4.27| ), are numerically integrated over one 
period of motion from the initial conditions, 

u{t = 0) = 1 
v{t = 0) = 

to obtain the resulting perturbation amplitudes at the end of one period of motion 
over the surface, u(T) and v{T). These quantities are compared with the general 
Floquet solution, ( [4.28| ), to obtain the equations, 

1 = a/i + bgi 

= af2 + bg2 . . 

u{T) = aX,h + bX29i ^^""^ 

v(T) = aAi/2 + bX2g2 

in which Ai = e'^^^ and A2 = e°^^"^ are the Floquet multipliers, and the Floquet 
functions are evaluated at t = (or t = T). The time symmetry relations, (|4.29|) , are 
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used to eliminate gi, (72, and A2 in ( |4.30| ), and the resulting equations are solved for 
Ai and ^, 



Ai,2 = u{T)±^uiTy-l 

The filament is linearly unstable to perturbations if and only if |tt(T)| > 1, with the 
growth parameter, 

a+ = ^ln(|S(T)| + v'ti(r)2-l) 
If u(T) > 1, the angle of this unstable mode, after each period, T, is 



9^ = arctan 



viT) 



and if u(T) < —1, 

-v{T) 



9a. = arctan 



u{Ty-l^ 



If |5(T)| < 1, the Floquet exponents are pure imaginary and the filament has neutral 
stability, with perturbations at this wavenumber oscillating, but not growing, about 
the filament center as it propagates over the surface. 



4.10 Stability results 

The wavy boundary significantly affects the evolution of the vortex filament. As a 
representative example, a plot of the growth parameter for a filament propagating 



over a boundary with height to length ratio -j- = .1 is shown in Figure |4]^. Previously 
stable wavenumbers become unstable, and the wavenumber and growth parameter 
of the maximally unstable mode increase. This indicates that the hairpin vortices 
evolving over a wavy boundary will be more tightly packed, having characteristic 
spanwise wavelengths less then lOOwu. The wavy boundary also has a significant 
effect on the growth angle of the hairpin vortices. The angle of the unstable modes. 



corresponding to the plot in Figure ^4.3| , is shown in Figure |4.4 The maximally 
unstable mode over this wavy boundary now grows at an angle of approximately 34°, 
a significant change to the 43° growth over a fiat plate. 

For very large values of the strong forcing is sufficient to induce new instability 
bands at previously stable wavenumbers, (Figure |475|). Although this new instability 
region has a smaller growth parameter than the primary instability, and will hence 
not significantly effect the evolution of a randomly perturbed vortex filament, it is 
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Figure 4.3: The growth parameter, 0"+, in units of 4, as a function of spanwise 
wavenumber, /cy*, for a filament over a wavy boundary of length L — 0.61mm, and 
height to length ratio, = 0.1. The maximum at ky^ ~ 0.072 is marked. The free 
stream velocity is U = 3—. The growth parameter for the fiat plate is shown for 
comparison, as a dotted curve. 
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Figure 4.4: The instability angle, 9^, as a function of spanwise wavenumber, ky^, for 
a filament over a wavy boundary of length L — 0.61mm, and height to length ratio, 
■| = 0.1. The maximally unstable wavenumber, kyi, ~ 0.072, is marked, indicating 
the angle of the maximally unstable mode, 9+ ~ 34°. The instability angle for the 
flat plate is shown dotted curve. 
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Figure 4.5: The growth parameter, cr+, in units of 4, as a function of spanwise 
wavenumber, A;^/^,, for a filament over a wavy boundary of length L = 0.61mm, and 
height to length ratio, 4 = 0.4. The free stream velocity is U = 3—. 



worth noting that the instability angle of the first of these new instability bands 
has the opposite sign of the previous instability, (Figure [4.6|) . A filament perturbed 
at this wavenumber grows in the previously stable direction, between —90° and 0°, 
illustrating the significant dynamical effect of the surface ridges. 

A contour plot of the growth parameter for a range of surface corrugation wave- 



lengths, L, is shown in Figure |4^, and a plot of the maximal growth angle is shown 
in Figure If the corrugations are small compared to the filament height, L <^vq, 
the surface is effectively smooth and does not affect the filament dynamics. Alter- 
natively, if the corrugation length is very large, L ^ vq, the parametric oscillations 
are too slow to affect the filament dynamics significantly. However, when the corru- 
gation length properly matches the filament height, L ~ 20vo, the dynamics of the 
filament are maximally affected, with an increase in maximally unstable wavenumber 
and a decrease in instability angle for the vortex filament. The surface corrugation 
wavelength of the dolphin, at L ~ .61mm, is well tuned to the corresponding height 
of vortex filament formation for a large range of swimming velocities. 

If the ridge length, L, and initial filament height in wall units, Vq, are taken as fixed 
properties of the boundary and developing flow, the instability may be calculated for 
a variety of ridge heights, ^, over a range of background velocities, U . This gives 
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Figure 4.6: The instability angle, 9+, as a function of spanwise wavenumber, ky*, for 
a filament over a wavy boundary of length L — 0.61mm, and height to length ratio, 
i = 0-4. 
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Figure 4.7: A contour plot of the growth parameter, -jj^, as a function of spanwise 
wavenumber, ky^, and corrugation wavelength, L, for a filament in background flow, 
U — 3—, over a wavy boundary of height to length ratio, 4 = 0.1. The corrugation 



wavelength for a dolphin, L ~ .61mm, is marked by a vertical line. The maximally 
unstable wavenumber, fc+y*, is plotted over the graph as a function of L. 
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Figure 4.8: A plot of the angle of the maximally unstable mode, 6*+, as a function 
of corrugation wavelength, L, for a filament in background flow, U — 3—, over a 
wavy boundary of height to length ratio, f — 0.1. The corrugation wavelength for a 
dolphin, L ~ .61mm, is marked by a vertical line. 
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Figure 4.9: The growth parameter, ^^p, as a function of spanwise wavenumber, 
ky^, and background velocity, U, for a filament at height Vq = Qwu propagating over 
a wavy boundary of height to length ratio = 0.2 and length L = .61mm. 



an indication of the effect of fixed length surface ridges for a range of swimming 
velocities. For delphin ridge spacing of L = .61mm and an initial filament height 
of Vq = 6wu, the instability and maximum instability angle is plotted, for ridges 
of height = .2, .1, .05, .025, over a large velocity range, (Figures |0| - |4.16| ). The 
maximum instability angle has a minimum at f/ ~ 5y - a velocity that is roughly 
independent of ridge height. The maximum instability angle at f/ ~ 5y ranges from 
27° for large ridges, = .2, to 40° for very small ridges, j = .025. 
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Figure 4.10: The angle of the maximally unstable mode, 9+, as a function of back- 
ground velocity: U, for a filament at height vq = 6wu propagating over a wavy 
boundary of height to length ratio j — 0.2 and length L — .61mm. 
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Figure 4.11: The growth parameter, -j^, as a function of spanwise wavenumber, ky^, 
and background velocity, U, for a filament at height vq = 6wu propagating over a 
wavy boundary of height to length ratio ^ = 0.1 and length L — .61mm. 
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Figure 4.12: The angle of the maximally unstable mode, as a function of back- 
ground velocity: U , for a filament at height vq = 6wu propagating over a wavy 
boundary of height to length ratio ^ = 0.1 and length L — .61mm. 



60 




2 4 e 8 10 12 14 



U (m/s) 

2 

Figure 4.13: The growth parameter, -jfj^, as a function of spanwise wavenumber, ky^, 
and background velocity, [/, for a filament at height vq = 6wu propagating over a 
wavy boundary of height to length ratio ^ = 0.05 and length L = .61mm. 



61 



45 



44- 
43- 
42- 



e+ 41 
40 



39 
38 
37 
36 



35 



10 



12 



14 



U (m/s) 



Figure 4.14: The angle of the maximally unstable mode, as a function of back- 
ground velocity: U , for a filament at height vq = 6wu propagating over a wavy 
boundary of height to length ratio ^ = 0.05 and length L — .61mm. 
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Figure 4.15: The growth parameter, -j^^, as a function of spanwise wavenumber, 
and background velocity, U, for a filament at height I'o = 6wu propagating over a 
wavy boundary of height to length ratio ^ = 0.025 and length L = .61mm. 
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Figure 4.16: The angle of the maximally unstable mode, as a function of back- 
ground velocity: U , for a filament at height I'o = &wu propagating over a wavy 
boundary of height to length ratio ^ = 0.025 and length L — .61mm. 
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Chapter 5 

The dolphins' secret 



By altering the dynamics of vortex filaments, dolphin skin reduces the rate of energy 
transport to turbulence in the evolving turbulent boundary layer. This results in 
a net reduction of the drag induced by skin friction. The effectiveness of this drag 
reduction technique is limited by the countering influence of increased skin friction 
due to the form drag of the surface ridges. Above a critical ridge height, j , the 
turbulent flow over the ridges separates, forming large eddies behind each ridge, and 
generates a large drag increase. The drag reduction at a given velocity will be optimal 
at ridge heights just below this critical value. 



5.1 The effect of surface corrugations 

The pondermotive forcing induced by the wavy boundary causes the vortex fllaments 
over the ridges to stretch at a maximally unstable angle, 6*+, that is significantly 
less then the maximally unstable angle for filaments over a flat boundary. This 
dynamical effect acts over a range of velocities centered about an optimum velocity 
determined by the wavelength of the surface corrugations. For a dolphin with surface 
ridges of wavelength L ~ .61mm, the optimal swimming speed for this dynamic 
effect is U ~ 5 — . Although most effective at this velocity, surface ridges of this 
wavelength produce a signiflcant reduction in the maximally unstable angle for a 
range of velocities between 2 — 15—, (Figures ^^[4.1(j|). This result is in excellent 



agreement with the reported swimming speed for this species of dolphin. 

The boundary layer vortex fllaments developing over such ridges grow into hairpin 
vortices inclined at the reduced angle. Packets of these hairpin vortices, now inclined 
at an angle lower than the 43° characteristic of hairpins over a flat surface, produce 
a signiflcantly lower Reynolds stress within the turbulent boundary layer. 
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5.2 Implications for turbulent energy transport 



The Reynolds stress contribution from hairpin vortices inclined at an angle 6^, (|2.8|) , 
is proportional to sin(26'+). Any deviation in hairpin inclination away from 45° pro- 
duces a reduction of Reynolds stress, with a corresponding reduction in the rate of 
turbulent energy transport. Since the dominant contribution to turbulent energy 
transport comes from the hairpin vortices, and the rate of energy transport is di- 
rectly proportional to the drag force on the surface, the percentage drag reduction 
due to the altered filament dynamics, can be roughly equated with the percentage 
reduction in Reynolds stress, 

sin (2 X 43°) - sin (29.) 
d. ~ 100- ^ ' ^ ^' 



sin (2 X 430) 

The drag reduction obtained using surface ridges is directly proportional to the degree 
to which the maximally unstable angle of the vortex filament may be lowered. This 
angle depends upon the swimming velocity and corrugation wavelength, as well as on 
the ridge height, ^. However, for ridge heights that are very large, the turbulent flow 
over the ridges will separate, and form large stationary eddies behind each ridge that 
effectively lower the ridge height and signiflcantly increase drag. 



5.3 Turbulent flow separation 

The drag reduction is limited by the maximal ridge height allowed before turbulent 
flow separation is induced. A linear model of turbulent boundary layer flow over a 
sinusoidal surface gives an estimate for the critical ridge height, ^ , at which sepa- 
ration occurs. This estimate is calculated, and confirmed by direct measurement, in 
papers by Zilker and Hanratty ||29|| , pO|] . For low friction velocities, Rei, = ^ < 630, 
and hence for low background velocities, the critical ridge height is calculated to be 
J ~ .033. However, the critical ridge height increases with velocity, because of the 
tendency of turbulent flow to inhibit separation. For example, at Re^^ ~ 6300, the 
critical ridge height becomes j ~ .1. Also, because the shape of cutaneous ridges 
over dolphin skin are not exactly sinusoidal, but have slightly steeper peaks and wider 
troughs, the critical ridge height for separation to occur over these ridges may vary 
from the value for perfectly sinusoidal ridges. It is likely that the shape and distribu- 
tion of the dolphins' cutaneous ridges is adapted to maximize the ridge height, and 
the dynamical effect on the vortex fllaments, without inducing flow separation. Given 
the structure of cutaneous muscle beneath the ridges, it is also very likely that the 
dolphins adjust the height of their cutaneous ridges to match the swimming velocity 
- keeping the ridge height slightly below the critical value in order to maximize drag 
reduction. 
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Figure 5.1: The estimated drag reduction due to altered filament dynamics, d-, 
plotted as a function of background velocity, U, for wavy boundaries of length L = 
.61mm and height to length ratio ^ = .1, = .06, and ^ = .033. 



5.4 Drag reduction 

The percent drag reduction due to the reduced inclination angle of the vortex instabil- 



ity for a variety of surface ridge heights is plotted in Figure |5Tl|. The drag reduction at 



the optimum swimming velocity of 5y varies from 2% for very small ridges of height 
J = .033 to 8% for ridges of height j = .1. Since it is likely that a dolphin increases 
the ridge height with swimming velocity, the percent drag reduction achieved is likely 
to increase with swimming velocity as a dolphin adjusts the ridge height to adapt to 
the higher critical ridge height. 



5.5 Application 

Man has long looked to the animal world for inspiration in efficient engineering de- 
sign. The discovery and analysis of cutaneous ridges in odontocetes presents another 
engineering feat of mother nature that may be emulated and put to human use. 

The drag reduction obtained through the use of surface corrugations will only be 
significant if properly matched to the selected application. The ridge spacing, L, must 
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be tuned for a specific velocity range, and the ridge height must not be so great as 
to induce separated fiow between the ridges. Although a dolphin may vary the ridge 
height to accommodate greater speeds, it is unlikely that an artificial implementation 
of variable ridge height would be cost effective to manufacture and employ. Therefore, 
cutaneous ridges may find practical application only when surfaces are expected to 
move through a fiuid at a steady velocity. 

The optimal ridge spacing for fiuid fiowing over a fiat plate is calculated to be 
L ~ 120j/*. This spacing depends only on the near boundary structure of the turbulent 
boundary layer, particularly the length of one wall unit, y*. The length of a wall unit 
over a smooth boundary may be determined experimentally by measuring the mean 
surface stress, tq, by observing the spanwise wavelength of the coherent structures 
present, ~ lOOy*, or, alternatively, by approximating a wall unit as ^ 27^. For 
example, an aircraft which regularly travels at a speed of 150— may wish to employ 
surface corrugations of wavelength 

/ 1.5xl0-5^\ 
L ~ 120 27 ~ .3mm 

V i50f ; 

The ridge height should be set just below the critical value, to avoid separation. 
A ridge height ratio of ^ ~ .03 will always be below the critical value, and produce a 
surface drag reduction of 2%. However, the critical value of the ridge height should be 
determined experimentally, so that the maximal ridge height may be employed and 
the greatest drag reduction achieved. It is possible that ridge heights of ratio greater 
than .1 may be allowed, producing an 8% or higher reduction in surface drag. 



5.5.1 Experiment 

Although several researchers have investigated the development of a turbulent bound- 
ary layer over a wavy wall, through experiment and numerical simulation |p9|] pO| ||3l| ||3^ 



none have focused on the parameter range relevant to the fiow of water over dolphin 
skin, or given attention to the development and evolution of coherent structures. 
However, it should be straightforward to extend these investigations to the relevant 
parameter range and thus determine the validity of the theoretical predictions. 

The optimal ridge height for drag reduction, the critical height at which separation 
occurs, should be determined experimentally. It would also be very usefull to deter- 
mine if a swimming dolphin does increase the cutaneous ridge height with swimming 
velocity, and to what degree; although it is not clear how such an experiment would 
be carried out. The dolphins have had sixty million years to optimize the reduction 
of drag through the control of vortex filament dynamics, and there is still a great deal 
to learn from these fascinating creatures. 
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